INTEGERS WITHOUT LARGE PRIME FACTORS: 
FROM RAMANUJAN TO DE BRUIJN 



PIETER MOREE 
In memoriam: Nicolaas Govert ('Dick') de Bruijn (1918-2012) 



Abstract. A small survey of work done on estimating the num- 
ber of integers without large prime factors up to around 1950 is 
provided. Around 1950 N.G. de Bruijn published results that dra- 
matically advanced the subject and started a new era in this topic. 

1. Introduction 

Let P{n) denote the largest prime divisor of n. We set -P(l) = L A 
number n is said to be y-friahl(Q if P{n) < y. We let S{x, y) denote the 
set of integers 1 < n < x such that P{n) < y. The cardinality of S{x, y) 
is denoted by ^{x, y). We write y = x^l'"" , that is m = logx/ logy. 

Fix M > 0. In 1930, Dickman proved that 



(1) hm ^^^'^ ^ =p{n\ 

with 

p(u) = p(N) - f ~ ^h v, (N <u<N+l, = 1, 2, 3, . . .), 
Jn V 

and p{u) = 1 for < M < 1 (see Figure 1). It is left to the reader to 
show that we have 



(2) p{u) 



1 for < n < 1; 

i p{u — t)dt for u> 1. 

The function p{u) in the literature is either called the Dickman function 
or the Dickman- de Bruijn function. 

In this note I will briefly discuss the work done on friable integers 
up to the papers of de Bruijn [71 18] that appeared around 1950 and 
dramatically advanced the subject. A lot of the early work was carried 
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out by number theorists from India. 

De Bruijn [7] improved on ([T]) by establishing a result that together 
with the best currently known estimate for the prime counting function 
(due to I.M. Vinogradov and Korobov in 1958) yields the following 
result. 

Theorem 1. The estimate 

(3) *(,,,)=,,(„){i+ca(M^)}, 

holds for 1 <u < log^^^'^y, that is, y > exp(log^'"^^^ x). 



De Bruijn's most important tool in his proof of this result is the 
Buchstah equation J9j, 

(4) ^!{x,y) = -^{x,z)- V 

^ p 

y<P<z 

where 1 < y < z < x. The Buchstab equation is easily proved on noting 
that the number of integers n < x with P{n) = p equals "^{x/pyp). 
Given a good estimate for \E'(x, for u < h, it allows one to obtain a 
good estimate for u < h + 1. 

De Bruijn [8] complemented Theorem [1] by an asymptotic estimate 
for p{u). That result has as a corollary that, for m > 3, 

log2_u - 1 , ^^^log2M^2^ 
logM 

which will suffice for our purposes. Note that ([5]) implies that, as 
U — 7- oo, 

formulas that suffice for most purposes and are easier to remember. 
For a more detailed description of this and other work of de Bruijn in 
analytic number theory, we refer to Moree |20] . 



(5) „„).e.p{-.{log..,og,.-l.!2|^-.0((i;|^)^)}}, 



2. Results on p{u) 

Note that p{u) > 0, for if not, then because of the continuity of p{u) 
there is a smallest zero uq > 1 and then substituting uq in ([2]) we easily 
arrive at a contradiction. Note that for m > 1 we have 

(6) p iu) = 

u 



^The reader not familiar with the Landau-Bachmann 0-notation we refer to 
wikipedia or any introductory text on analytic number theory, e.g., Tenenbaum 
[35] , Instead of log log a; we sometimes write logaX, instead of {\ogx)^, log'^a;. 
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It follows that p{u) = 1 — logu for 1 < u < 2. For 2 < u < 3, 
p{u) can be expressed in terms of the dilogarithm. However, with 
increasing u one has to resort to estimating p{u) or finding a numerical 
approximation. 

Since p{u) > we see from that p{u) is strictly decreasing for 
u > 1. From this and ([2]) we then find that up{u) < p{u — 1), which on 
using induction leads to p{u) < l/[u]\ for u > 0. It follows that p{u) 
quickly tends to zero as u tends to infinity. 

Ramaswami [22] proved that 

c 

for a suitable constant C, with F the Gamma function. By Stirling's 
formula we have log F(m) ~ u \ogu and hence the latter inequality is for 
u large enough improved on by the following inequality due to Buchstab 

(7) piu) > exp I — n| log u + logo u + 6 ^ 1 1 , iu>6). 

I I logM J J 

Note that on its turn de Bruijn's result considerably improves on 
the latter inequahty. 

3. S. Ramanujan (1887-1920) and the friables 

In his first letter (January 16th, 1913) to Hardy (see, e.g. |3j), one of 
the most famous letters in all of mathematics, Ramanujan claims that 

(8) v,(n,3) = l-l°S(2n)log(3n) 



2 log 2 log 3 

The formula is of course intended as an approximation, and there is no 
evidence to show how accurate Ramanujan supposed it to be. Hardy 
[T7t pp. 69-81] in his lectures on Ramanujan's work gave an account 
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of an interesting analysis that can be made to hang upon the above 
assertion. I return to this result in the section on the \E'(a;, y) work of 
Pillai. 

In the so-called Lost Notebook [27] we find at the second half of page 
337: 

is the no. of nos of the form 

2'^2 .3a3 .5a5...^ap ^ ^ 

not exceeding x. 



^ < e < 1, 



2 



A, 







4 



(x) ~ x| 1 



1 <e< i (hix) -xU - I + 



Ai Ai Ao J 

^ rfAo n dXi f^-^' dXo 



Ao ie Ai Jx^ Xq 

^ t/Aa dXi f^~^' dXc 



A2 J\n Ai A Ao 



1 , , 1 , . ^ I , dXo dXi dXo 



Ao Je Ai 
3 rfAa /"^ rfAi /"^-^^ dXo 



+ 



A2 JAa -^1 ^Ai Ao 

3 rfAs /"^ dX2 dXi f^'^^ dXo 



A3 A3 ^2 A2 -^1 Ai Ao 



and so on. 



In the book by Andrews and Berndt [Tj §8.2] it is shown that Ra- 
manujan's assertion is equivalent with ([1]) with 

fc=0 

where 

/■ dti dtk 
hiu) = — . 

J H,-tk>i ti tk 

ti+... + tf.<u 

This is one of many examples where Ramanujan reached with his hand 
from his grave to snatch a theorem, in this case from Dickman who 
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was at least 10 years later than Ramanujan, cf. Berndt [2]. Chowla 
and Vijayaragahavan [12] seemed to have been the first to rigorously 
prove ([1]) with p{u) expressed as a sum of iterated integrals (cf. the 
section on Buchstab). The asymptotic behaviour of the integrals Ik{u) 
has been studied by Soundararajan [32] . 

Ramanujan's claim reminds me of the following result of Chamayou 
[To]: If xi,X2,X3, ■ ■ ■ are independent random variables uniformly dis- 
tributed in (0, 1), and Un = Xi + X1X2 + . . . + X1X2 ■ ■ ■ Xn, then m„ 
converges in probability to a limit Uoo and Uoo has a probability distri- 
bution with density function p{t)e~'^, where 7 denotes Euler's constant. 

4. I.M. Vinogradov (1891-1983) and the friables 

The first to have an application for \l'(x, estimates seems to have 
been Ivan Matveyevich Vinogradov [SH] in 1927. Let > 2 be a pre- 
scribed integer and p = l(mod k) a prime. The k-th powers in (Z/pZ)* 
form a subgroup of order (p — l)/k and so the existence follows of 
gi{p, k), the least k-th power non-residue modulo a prime p. Suppose 
that y < gi{p, k), then S{x,y) consists of k-th power residues only. It 
follows that 

^{x,y) < 4f{n < X : n = a'^(mod p) for some a}. 

The idea is now to use good estimates for the quantities on both side 
of the inequality sign in order to deduce an upper bound for gi{p, k). 

Vinogradov [36] showed that '^{x^x^^'^) > 5{u)x for x > 1, m > 0, 
where 5{u) depends only on u and is positive. He applied this to show 
that if m > 8, /c > m™', and p = l(mod k) is sufficiently large, then 

(9) g^{p,k)<p'/^. 

See Norton [21] for a historical account of the problem of determining 
gi (p, k) and original results. 

5. K. DiCKMAN (1861-1947) and the friables 

Karl Dickman was active in the Swedish insurance business in the end 
of the 19th century and the beginning of the 20th century. Probably, 
he studied mathematics in the 1880's at Stockholm University, where 
the legendary Mittag-Leffier was professoil^. 

As already mentioned Dickman proved ([T]) and in the same papei@ 



I have this information from Lars Hoist. 

^Several sources falsely claim that Dickman wrote only one mathematical paper. 
He also wrote |14) . 
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gave an heuristic argument to the effect that 




du. 



Denote the integral above by A. Diclcman argued that A ~ 0.62433. 
Mitchell [19] in 1968 computed that A = 0.62432998854 .... The in- 
terpretation of Dickman's heuristic is that for an average integer with 
m digits, its greatest prime factor has about Am digits. The constant 
A is now known as the Golomb-Dickman constant, as it arose indepen- 
dently in research of Golomb and others involving the largest cycle in 
a random permutation. 

De Bruijn in 1951 was the first to prove fllOp . He did this using 
his A(a;, |/)-function, an approximation of \E'(x,y), that he introduced 
in the same paper. 

6. S.S. PiLLAI (1901-1950) AND THE FRIABLES 

Subbayya Sivasankaranarayana Pillai (1901-1950) was a number the- 
orist who worked on problems in classical number theory (Diophantine 
equations, Waring's problem, etc.). Indeed, he clearly was very much 
inspired by the work of Ramanujan. He tragically died in a plane crash 
near Cairo while on his way to the International Congress of Mathe- 
maticians (ICM) 1950, which was held at Harvard University. 

Pillai wrote two manuscripts on friable integers, [23l [23] , of which 
[23j was accepted for publication in the Journal of the London Math- 
ematical Society, but did not appear in print. Also [23| was never 
published in a journal. 

In [23], see also [2S1 PP- 481-483], Pillai investigates 'if{x,y) for y 
fixed. Let pi,p2, . . . ,Pk denote all the different primes < y. Notice 
that \l/(a;, y) equals the cardinality of the set 



Thus \l/(x,?/) equals the number of lattice points in a /c- dimensional 
tetrahedron with sides of length loga;/log2, . . . , logx/logpfc. This 
tetrahedron has volume 



k 



{(ei, . . . , Cfc) G Z'' : Ci > 0, ^ logpi < x}. 




Pillai shows that 
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If pi, . . . , Pk are positive real numbers and pi/p2 is irrational, then the 
same estimate with logpj replaced by pi holds for 

k 

{(ei, . . . , Cfc) G Z'' : > 0, eipi < x}. 

i=l 

This was proved by Specht [33] (after whom the Specht modules are 
named), see also Beukers [1]. A much sharper result than that of 
Pillai/Specht was obtained in 1969 by Ennola [15] (see also Norton [211 
pp. 24-26]). In this result Bernoulli numbers make their appearance. 
Note that Pillai's result implies that 

/I IX ,T,/ ON 1 log(2x) log(3x) 

(11) ^^^'^) = 2 log21og3 

and that the estimate 

2) = o/"oi^ Q + ^^^^^ ^) 
2 log 2 log 3 

is false. Thus Ramanujan's estimate (|8]) is more precise than the trivial 
estimate log^x/(21og21og3). Hardy [LT^ §5.13] showed that the error 
term o(logx) in ( TTTl) can be replaced by o(logx/ log2 x). In the proof 
of this he uses a result of Pillai [22], see also ^3 PP- 53-61], saying 
that given < 5 < 1, one has |2^ — 3^| > 2*^^^"^^^ for all integers x and 
y with X > xq{6) sufficiently large. 

In [2^, see also pp. 515-517], Pillai claims that, for u > 6, 
B/u < p{u) < A/u^ with Q < B < A constants. He proves this result 
by induction assuming a certain estimate for p(6) holds. However, this 
estimate for p(6) does not hold. Indeed, the claim contradicts ([5]) and 
is false. 

Since Pillai reported on his work on the friables at conferences in In- 
dia and stated open problems there, his influence on the early develop- 
ment of the topic was considerable. E.g., one of the questions he raised 
was whether ^!{x^x^/^) = 0{x^^^) uniformly for u < (logx)/ log2. 
This question was answered in the affirmative by Ramaswami 



7. R.A. Rankin (1915-2001) and the friables 

In his work on the gaps between consecutive primes Robert Alexan- 
der Rankin [3T] in 1938 introduced a simple idea to estimate \l/(x, y) 
which turns out to be remarkably effective and can be used in simi- 
lar situations. This idea is now called 'Rankin's method' or 'Rankin's 
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trick'. Starting point is the observation that for any a > 

(12) ^ix,y)< Yl C-r<^'' E ^.=^''c{o^yi 

n£S{x,y) P{n)<y 

where 

p<y 

is the partial Euler product up to y for the Riemann zeta function 
Recall that, for 3fJs > 1, 



- V 

n=l p 



By making an appropriate choice for a and estimating C{cr,y) using 
analytic prime number theory, a good upper bound for \E'(x, can be 
found. E.g., the choice a = 1 — 1/(2 logy) leads to 

aa,y) « exp { E ^} ^ { E E ^) } « 1°^^' 

p<y p<y p<y 

which gives rise to 

(13) "^{x, y) < xe""/^ log y. 



8. A. A. BuKHSHTAB (1905-1990) and the friables 

Aleksandr Adol'fovich Bukhshtab (Buchstab in the German spelling) 
most important contribution is the equation now named after him. 
A generalization of it plays an important role in sieve theory. Buchstab 
[H] in 1949 proved ([T]) and gave both Dickman's differential-difference 
equation as well as the result 
(14) 

p{u) = 1 + E(-i)" r T"' r~' ■ ■ ■ T""' dudu.,...du ^ 

Jn Jn~l Jn-2 J I ^1*^2 " " " i^n 

for < M < A^ + 1 and > 1 an integer, simplifying Chowla and Vija- 
yaragahavan's expression (they erroneously omitted the term n = N). 
Further, Buchstab established inequality ([7]) and applied his results 
to show that the exponent in Vinogradov's result can be roughly 
divided by two. 
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9. V. Ramaswami and the friables 
V. Ramaswamll [28] showed that 

X , 



^^/{x,x^/'^) = p{u)x + OuU 



logx 

for X > 1, 1 < u < U, and remarked that the error term is best possible. 
He sharpened this result in [29] and showed there that, for u > 2, 

(15) =p(u)x + (T(n)-^ + 0( ^ ^ 



logx log^/^x 

with (t{u) defined similarly to p{u). Indeed, it turns out that 

a{u) = (1 -7)p(u- 1), 

but this was not noticed by Ramaswami. In [30] Ramaswami gener- 
alized his results to Bi{m,x,y) which counts the number of integers 
n < X with P{n) < y and n = /(mod m) [^. Norton [2ll pp. 12-13] 
points out some deficits of this paper and gives a reproof [211 §4] of 
Ramaswami's result on i?;(m, x, x^/") generalizing (fT5l) . 

From de Bruijn's paper [TJ Eqs. (5.3), (4.6)] one easily derives the 
following generalization of Ramaswami's result^: 

Theorem 2. Let m > 0, x > 1, and suppose m + 1 < u < y/log x . 
Then 

'^{x,y) =xVa^ ^ ["^ +Qm( , ^+1 ), 
log y Mog™+'i/' 

with the r-th derivative of p{u) and oq, ai, . . . are the coefficients 

in the power series expansion 

Z 

C(l + z) = ao + aiz + a2Z^ + . . . , 

\ ~\~ z 

with \z\ < 1. 

It is well-known (see, e.g., Briggs and Chowla [5]) that around s = 1 
the Riemann zeta function has the Laurent series expansion 



A;=0 

with 7fc the /c-th Stieltjes constant and with 70 = 7 Euler's constant. 
Using this we find that = 1 and ai =7 — 1. Thus Theorem [2] yields 



^He worked at Andhra University until his death in 1961. I will be grateful for 
further bibliographical information. 

^Buchstab [9j was the first to investigate Bi{m,x,y). 

'''The notation Om indicates that the implied constant might depend on to. 
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( flSll with a(u) = (1 — 'y)p{u — 1) for the range 2 < u < ^/\ogx. For 
u > y/logx the estimate f|T5l) in view of ([5]) reduces to 

which easily follows from f|T3l) . 

10. S. Chowla (1907-1995) and the friables 

The two most prominent number theorists in the period following 
Ramanujan were S.S. Pillai and Sarvadaman Chowla. They kept in 
contact through an intense correspondence [M]. Chowla in his long 
career published hunderds of reseach papers. 

Chowla and Vijayaragahavan [12] expressed p{u) as an iterated in- 
tegral and gave a formula akin to (IT^ . De Bruijn [6] established some 
results implying that ^{x,log^ x) = 0{x^~^^^~^^) for h > 2. An easier 
reproof of the latter result was given by Chowla and Briggs |llj . 

11. Summary 

It seems that the first person to look at friable integers was Ra- 
manujan, starting with his first letter to Hardy (1913), also Ramanujan 
seems to have been the first person to arrive at the Dickman-de Bruijn 
function p{u). Pillai generalized some of Ramanujan's work and spoke 
about it on conferences in India, which likely induced a small group 
of Indian number theorists to work on friable integers. Elsewhere in 
the same period (1930-1950) only incidental work was done on the 
topic. Around 1950 N.G. de Bruijn published his ground-breaking pa- 
pers [II l8j. Soon afterward the Indian number theorists stopped pub- 
lishing on friable integers. 

It should also be said that the work on friable integers up to 1950 
seems to contain more mistakes than more recent work. Norton |2l] 
points out and corrects many of these mistakes. 

Further reading. As a first introduction to friable numbers I can 
highly recommend Granville's 2008 survey [12]. It has a special em- 
phasis on friable numbers and their role in algorithms in computational 
number theory. Mathematically more demanding is the 1993 survey by 
Hildebrand and Tenenbaum [TS]. Chapter III. 5 in Tenenbaum's book 
[35] deals with p{u) and approximations to \l/(x, y) by the saddle point 
method. 

Acknowledgement. I thank R. Thangadurai for helpful correspondence 
on S.S. Pillai and the friables and providing me with a PDF file of 
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Pillai's collected works. B.C. Berndt kindly sent me a copy of K.K. 
Norton provided helpful comments on an earlier version. His research 
monograph [21], which is the most extensive source available on the 
early history of friable integer counting, was quite helpful to me. In 
plj . by the way, new results (at the time) on gi{p, k) and "^mix, y), the 
number of ^/-friable integers 1 < n < x coprime with m are established. 
Figure 1 was kindly created for me by Jon Sorenson and Alex Weisse 
(head of the MPIM computer group). 
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